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Abstract 

Consolidation of colloidal gels results in enormous stresses that are 
usually released in the formation of undesirable cracks. The capacity of 
a gel network to crack during drying depends on the existence and sig- 
nificance of a pressure gradient in the pore liquid; in addition it depends 
on the way the gel relaxes the resulting drying stresses. In this paper the 
effect of a binary mixture of solvents saturating the gel network on the 
crack patterns formation is investigated. Indeed, incorporation of a small 
quantity of non-volatile cosolvent, i.e. glycerol, inhibits drying-induced 
cracks; moreover addition of a concentration greater than 10% to a col- 
loidal dispersion leads to a crack free coating in room conditions. Mass 
variation with time reveals that both evaporation rate and cracking time 
are not affected by glycerol, in the range of added glycerol contents stud- 
ied. In addition measurements of mechanical properties show that the 
elastic modulus is reduced with glycerol content. The decrease of the 
number of cracks with the glycerol content is related to the flattening of 
the pressure gradient in the pore liquid. The mechanism is shown to be 
due to the combination of two processes: flow driven by the pressure gra- 
dient and diffusion mechanisms in accordance with Scherer work (1989). 

1 Introduction 

Most coatings are made by depositing a volatile liquid that contains dispersed 
colloidal particles on a surface. The liquid is then evaporated until a dry film is 
obtained. Common examples in industrial fields are paints, protective coatings 
(anti-corrosion), ceramic membranes [3]. Coatings are required to be homoge- 
neous and smooth. Particularly they have to be free of heterogeneities possibly 
leading to cracking due to high stresses developing during the consolidation. 
In a liquid film, evaporation of the volatile solvent causes the solid particles 
to be confined into a smaller volume, until they come into direct contact with 
each other, and form a solid porous network, so called gel !4j. Then, further 
evaporation causes the liquid/air interface to curve into a set of menisci that 
join the particles. As a consequence the capillary forces exert extremely strong 
compressive stress on the particle network. When the drying stress development 
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in the gel exceeds the strength of the material, undesirable effects occur such as 
the formation of various modes of cracking [Tl [2"Tl |3"T1 I3"2"] . As a consequence the 
stress in the gel is relaxed. Since numerous studies on cracks induced by dry- 
ing of colloidal suspensions have been investigated for the last two decades, the 
stress development and the organization of colloidal particles in a drying film 
remain unclear. One of the main questions that still stands: how can the stress 
be relaxed before cracks can form? In addition, the development of stresses, 
hence of cracking notably depends on the type of the particles (size[T2], chem- 
ical compositional; 26 ) and mechanical properties of the gel formed f[2l \TE\). 
the substrate, the drying conditions ftemperature|15j . relative humidity [11], air 
velocity[17]) and the film thickness [2l [TBI 17] . In particular the capacity of a gel 
network to crack during drying depends on the existence and significance of a 
pressure gradient in the pore liquid and the way the gel relaxes the induced 
drying stresses. 

In the present work, we focus on the effect of cosolvent added to the col- 
loidal dispersion on crack patterns induced by drying. Such additive reveals 
crack- free drying when added to ceramic foams [T4"]. Glycerol has been chosen 
as a cosolvent for interesting properties of high miscibility in water, low volatil- 
ity, and because it exhibits a highly different intrinsic diffusion coefficient from 
that of water. Also glycerol is an element composing many essential fine arts 
products [16]: due to its transparency properties and because it is used as a 
softener, avoiding creases and cracks. Crack patterns have been quantified as 
a function of the added glycerol content. Drying kinetics have been investi- 
gated by mass variation with time measurements. Also mechanical properties 
of colloidal gels are measured by indentation testing as a function of the added 
cosolvent content. Experimental results can be explained by the combination of 
two processes occuring during the drying process: flow driven by the pressure 
gradient accordingly with Darcy law and diffusion mechanism accordingly with 
Fick law. 

2 Experimental 

2.1 Materials and samples preparation 

We used colloidal silica dispersion, Ludox HS-40 purchased from Sigma- Aldrich. 
The radius of the particles is 9 ± 2 nm, the density 2.36 x 10 3 kg/m 3 . The 
mass fraction of the initial dispersion was estimated to <p m — 40.5% by a dry 
extract. Thus, we could deduce the particle volume fraction for this solution 
to be 4> v ~ 22%. Particles are polydisperse enough not to crystallize. In the 
absence of evaporation, the suspension is stable due to the competition between 
the attractive Van der Waals forces and the electrostatic repulsion [20] between 
particles (pH = 9.8). 

Glycerol (purity: 99.0%, from Merck) was added to the aqueous dispersion 
as a cosolvent. This solvent exhibits high miscibility properties in water and low 
vapor pressure (2.2 x 10~ 4 mbar at 25°C)[pJ compared to water (31.7 mbar). 
Thus, we assume glycerol as a non- volatile compound. The surface tension 
of pure glycerol is 62.5mN/m at 20° C [55]. Glycerol/water weight ratio was 
prepared using the following process: The resulting solution was filtered to 
0.5 /im to eliminate eventual aggregates; the mass fraction is estimated by a 
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Figure 1: (a) Experimental set-up (side view): during solvent loss the layer 
consolidates, (b) Sketch of the solid network saturated by solvent: curvature 
of the solvent/air menisci occurs at the evaporation surface during solvent loss 
and exerts an enormous stress on this network. 

dry extract at 140° C. In order to keep the particle mass fraction constant for 
each sample, a given quantity of water is evaporated, then replaced by an equal 
mass of glycerol; the evaporation of water is achieved by stirring the aqueous 
dispersion at 35°C in an Erlenmeyer. Finally, we checked the final pH of the 
solution. 

We note k the ratio between the mass of glycerol added and the mass of the 
resulting dispersion. Six solutions have been investigated: k — 1.1, 2.1, 3.1, 5.1, 
7.2 and 8.9 % wt. This leads to glycerol fraction in solvent denoted by Xi (the 
subscript stands for initial). In the range of glycerol content added, the solvent 
viscosity increases until 1.3 mPa.s. 

2.2 Drying geometry 

Experiments are performed on films of fixed thickness. Also a small amount 
of the solution (initial weight rrii = 0.52 ± 0.02 g) is deposited inside a circular 
container (diameter =15 mm) whose bottom is a glass microscope slide carefully 
cleaned with ethanol and the lateral wall is a rigid rubber ring (figure [1] (a)). 
The contact line of the solution is quenched at the upper edge of the circular 
wall and remains there during the whole experiment. In this geometry the 
layer dries at room temperature (T = 23° C) under controlled relative humidity 
(Rh ~ 50 ± 2%). In this way, the desiccation takes place in the absence of 
convection in the vapor so that the evaporation is limited by diffusion of the 
solvent into the air. The drying kinetics is obtained using a scale (Sartorius) 
with an accuracy of 0.01 mg. When the gelled phase is formed, the layer exhibits 
an approximately constant thickness in a central region (figure [1] (a)). In this 
region, covering about 70% of the total surface area, the evaporation is assumed 
to be uniform. A typical value of the thickness far from the border is h ~ 0.7 
mm. The colloidal gel is transparent allowing us to observe easily the dynamics 
of crack patterns in the layer. Also the crack patterns formation is recorded 
using an interval timer with an AVT's Marlin camera positioned on the top of 
the sample. 
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2.3 Measurements of the mechanical properties of the gel 
phase 

To characterize the mechanical properties of the gelled layers relevant elastic 
modulus, E p , has been performed by indentation testing. This quantity has 
been determined using a CSM Instruments Micro Indentation Tester (MHT) 
with a four-sided pyramidal Vickers indenter. Indentation testing is a local 
method of mechanical investigation. The indenter, initially in contact with the 
surface of the gel, is driven in the material until a maximal 70 mN load with 
a loading speed 70 mN/min. The maximal force is held during 10 s. Then, 
the load is decreased until zero with the same speed. Figure [SJi shows a typical 
load/unload curve. The standard way of estimating the elastic modulus from 
the indentation load-displacement curve uses the initial slope of the unloading 
curve. Based on Sneddon work |30j . the following formula has been derived for 
the elastic modulus, E p , to the unloading slope [241 [25] : 

F - J* S 

In this equation, S is the slope of the unloading curve at the start of un- 
loading, and A is the projected area of contact between the indenter and the 
material at that point (inset in figure [S^,). Note that the contact between the 
indenter and the gel is always chosen in a location where the layer adheres on 
the substrate. In the case of successive measurements of the elastic modulus 
with time, indentation testing was performed in the same adhering region to 
obtain comparable values as effectively as possible. Particularly, the difficulties 
of these measurements lie in gel films without glycerol addition exhibiting the 
smallest size of adhering regions. Indeed it has been observed that glycerol 
reduces the capacity of the gel to delaminate. In addition indentation testing 
applied to drying gelled layers has to be carefully analyzed since the system is 
not homogeneous in the thickness and evolves with time. Also measurements of 
elastic moduli in gels by this way do not give absolute values but provide on the 
one hand a process of comparing mechanical properties of different materials, 
on the other hand a characterization of the time evolution of the mechanical 
properties of the system. Therefore measurements for the elastic modulus have 
been reproduced several times and performed on 2 or 3 samples; bars in figure 
[SJd take into account the minimum and maximum measurements obtained. 



3 Results 
3.1 Drying 

Once the layer is deposited in the circular trough, evaporation of the water 
takes place. The drying process of a colloidal dispersion is usually separated 
into two stages, that can be roughly distinguished by measuring mass, m, and 
drying rate, variations with time ( figure l2"j) |29j . In the first regime (con- 
stant rate period) the evaporation rate is constant and evaporation is mainly 
controlled by the external conditions of relative humidity and temperature in 
the surroundings. During this stage, water removal concentrates the colloidal 
particles into a closed packed array: a porous network saturated with solvent. 
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Figure 2: Drying curves. Drying rate, for silica films for different additional 
glycerol concentrations, n. The drying rate is calculated from the mass variation 

Ft 1 

with time (upper inset). The variation 4^ with time evidences the transient 
time close to 32000s. 

Due to geometry of the planar film, evaporation at the gel-air interface and 
adhesion at the gel-substrate interface result in the development of high stresses 
in the gel (figure QTb)). Capillary pressure that is caused by the liquid menisci 
formed between the particles at the top of the packed region is responsible for 
the shrinkage of the porous network. The shrinkage of the gel at the gel-air 
interface is frustrated by the adhesion to the substrate. As a consequence ten- 
sile stresses progressively build up in the layer. Then a non-linear drying rate 
follows: the falling rate period, when liquid/air interfaces recede into the porous 
medium and when the drying process is limited by flow of liquid to the top layer 
through the porosity. The transient between the constant rate period and the 
falling rate period is highlighted in the lower inset in Figure [5] it shows the 
drying rate derivative ^2™ . In our experimental conditions we observe that the 
transient always occurs at ~ 32000s in the studied k range. Consequently the 
drying rate decreases as clearly shown in figure [2] Similar behaviors occur when 
glycerol content is added to the aqueous dispersion (figure (2). 

3.2 Cracking 

As the tensile stress reaches a threshold value, cracks appear in the film and 
invade the plane of the layer. First cracks appear at t £ [24000, 27000]s after 
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Figure 3: Mass ratio of samples for different glycerol concentration n. The 
inset shows the final concentration of glycerol in the solvent wetting the porous 
medium as a function of k. The solid line is a linear adjustment (slope= 11.3 ± 
0.3). 



deposition of the layer on the substrate, without any correlation with the glycerol 
concentration nor the initial mass disparity. Moreover the cracks invade the 
layer during only a few minutes. Thus the crack network develops during the 
initial stage of the drying process, in constant rate period [3J and does not 
evolve afterwards. Plot in inset in figure [3] shows the ratio, m^/m^, of the 
final (at t — 60000s) and the initial mass of material on the glass plate (raw 
data). Futhermore, we know that silica particles contributes to 40% of the mass. 
Thus, the excess of mass is due to the solvent. As previously mentioned only 
water could evaporate and glycerol remains in the porous medium. Also we can 
calculate the mass fraction of glycerol \f composing the solvent when the mass 
variations with time do not evolve anymore (t ~ 60000s): 

_ mf/m t - 0.4- k 
Xf ~ nif/m, - 0.4 U 

We found that the higher the initial fraction k is, the higher Xf reaches (see 
inset of figure[3|), leading to a saturation of glycerol in the gel for n > 8.9 ±0.3% 
in accordance with the fit. 

Crack patterns are strongly modified by the presence of a glycerol within 
the solvent. Indeed, in the case of a silica gel without additional glycerol, the 
successive formation of cracks results in the typical pattern shown in figure 
[4^. This pattern is usually continued by further crack generations: particularly 
delamination process takes place when the gel detaches from the substrate [25] . 

Moreover final crack patterns are shown in figure |3Ji,b,c,d for increasing glyc- 
erol concentrations n up to 10.2% in our experimental conditions. For each pat- 
tern, cracks still interconnect leading to complete network: no broken network 
is observable suggesting that if crack nucleation is less favourable in presence 
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Figure 4: Crack patterns after 16 hours in gelled layers when adding a second 
solvent to the system, (a) Without any additional glycerol, (b) k = 5.1%. 
(c) K — 8.9%. (d) K — 10.2%: the layer is crack-free at the mesoscopic scale, 
(e) Crack spacing plotted as a function of the concentration in glycerol (drying 
conditions Rh — 50% and T = 23°C; the initial weight deposited in each 
trough is the same), (f) Propagating crack tip in a gel with additional glycerol 
k = 8.9%. 

of additional glycerol, crack propagation is still not inhibited. The number of 
cracks generated diminishes when the quantity of additional glycerol increases; 
this results in increasing the crack spacing as shown in figure^- FigureHf shows 
a magnification of a typical crack tip propagating in a gel with glycerol content. 
The sharpness of this crack tip, also observed for every k < 9.5%, suggests that 
the gel still exhibits brittleness properties. In addition buckle-driven delamina- 
tion process, that usually takes place after the formation of a crack network, is 
inhibited by increasing additional glycerol concentration in the gel. 

Above a threshold concentration of additional glycerol, the layer remains 
crack free (figure HJl) : no crack forms because the gel could not reach the com- 
paction (or the threshold stress) for which cracks propagate. It results in a flat 
gel film saturated in glycerol. 

3.3 Gel consolidation 

The macroscopic elastic response of the gel is characterized by indentation test- 
ing for various binary mixtures of solvents (figure 0). Based on the procedure 
depicted in figure [5^,, elastic moduli have been measured for gels of various glyc- 
erol contents. Each measurement has been obtained ~ 8 hours after the cracks 
formation since the gel evolves slowly after this duration (figure [5]) . Measure- 
ments in figure \5jp reveal that an increase of the quantity of glycerol in the gel 
results in a decrease of its elastic modulus. Since glycerol remains inside the 
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Figure 5: Macroscopic elastic response of the gel obtained by indentation test- 
ing, (a) Loading/unloading curve: force applied on the tip versus indcntcr 
displacement (pm). Insets: variation of the applied force versus time; image 
of indentation print at the surface of the gel after a loading/unloading cycle 
observed using optical microscopy (maximum load: 70mN, bar = 25/im): the 
dashed square limit the projected area of contact, A. (b) Variation of the elas- 
tic modulus, E p , for seven glycerol-Ludpx solutions, (c) Variation of the elastic 
modulus, Ep, as a function of time forXudox HS-40 without glycerol addition 
(k = 0%) and for two different concentrations of glycerol (k = 3.1%, 8.9%). 
Measurements start after crack propagation, t > t m . 



porous medium, it could act as a wetting agent lowering the material stiffness. 

The elastic modulus of silica gels has been measured as a function of time 
during the drying process; measurements started after crack propagation, t > 
t m . These measurements are shown for the silica gel without any glycerol ad- 
dition and for two glycerol contents in figure [5p. In any cases, gel films clearly 
stiffens as they consolidate with time; for the silica gel without any glycerol 
addition, this process is more rapid than in presence of glycerol content. Two 
factors could contribute to increase the modulus as a function of time: stiffening 
of the solid phase of the gel and reduction in the porosity [37] . Also the elastic 
modulus increases faster for gels having lower concentrations of glycerol. 

4 Discussion 

The capacity of a gel to crack during drying depends on the existence and sig- 
nificance of a pressure gradient in the pore liquid. Indeed the spatial variation 
in pressure causes a variation in contraction of the gel network |28) . This mod- 
ifies the way the drying stresses are relaxed and consequently the crack are 
formed. In the following the spatial variation in pressure in the liquid pore is 
characterized in gels saturated by binary mixtures of solvents during the first 
drying regime since cracks form during this stage. Particularly we have cho- 
sen to compare the pressure distribution for different glycerol contents at time 
t m = 25500s. At this moment, close to the typical cracking time of the gelled 
layers, gels can be considered as non-shrinking, that is the gel network does not 
change during time|28); pores contain the binary mixture of solvents and the 
liquid-air interfaces are already curved. 

In drying gels, liquid usually moves through the network in response to a 
pressure gradient of liquid, ^j, accordingly to the Darcy law[4]. Also, in the 
case of a gel saturated with a single volatile solvent, the flux of evaporating 
liquid at the gel-air interface comes entirely from flow of liquid. It expresses as: 



where k is the network permeability of the gel, and r\ is the dynamic viscosity 
of the solvent in the pore; note that the pressure of the liquid pore is regarded 
as positive. 

Let us now consider the incorporation of a soluble and non- volatile cosolvent, 
i.e. glycerol, to the aqueous dispersion. The drying process leads to a gel 
phase saturated by a binary mixture of water and glycerol. This results in a 
combination of the pressure flow component © and a diffusion mechanising, 28, 
3]. In that way, the diffusive flux for each component, i, expresses in accordance 
with the Fick law: 



where Di is the intrinsic diffusion coefficient refers in terms of the rate of 
transfer of component i across a fixed section. Note that P and Ci are time 
dependent. 

The total fluid flux in the porous network combines the flow driven by the 
pressure gradient @ and the diffusive flux ©. Also, the relative change in 
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volume of liquid over time, Vl/Vl, is due to the change in both local flow, 
d z Joarcy, and local diffusion flux, d z J l Fick . As a result the following conservation 
equation can be written [2%]: 

(1 - P)^r = -d z J D arcy ~ Kd z {J Flck + J Flck ) (4) 

where p is the relative density of the gel, K denotes the porosity of the gel, 
J Fick is the diffusive flux of component i, stating i = 1 for water and i = 2 for 
glycerol. 

Since the change in the volume of liquid has to equal the change in pore 
volume, the following relation can be written in accordance with reference [28]: 

(i - PKr = -if O + (i-P)<P>) (5) 



V L E. 



where (3 is a function of the Poisson ratio of the gel and P denotes the 
time derivative of the pressure. Also the flux from the gel surface equals the 
rate of change in volume of the gel; this expresses as < P >= E V ^P- for 
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layer of thickness h accordingly to reference [28]. Here the evaporation rate, Ve, 
determines the rate of volume loss by water evaporation. 

By equating the right sides of expressions 2] and [5] and taking into account 
relations [5] and [31 it comes: 

P dP , m V E kd 2 P r , .„ „ , d 2 d 
W p -8i + ^-^ = -^ + K ^- D ^^ (6) 

Equation [6] reveals that the diffusion mechanism can become significant only 
when the intrinsic diffusion coefficients of cosolvents are sufficiently different 
from each other. Note that equation [5] is only valid in the case of a solvent 
made of two components: in the case of a pure solvent, the diffusive term 
(second term of the right side) disappears and the flow is only governed by the 
Darcy law. Applying equation [B] to our system allows us to plot the pressure 
in a pore liquid of the gel for different glycerol contents. In that purpose the 
concentration of component 1 (water), C\, has to be incorporated in equation^ 
C\ satisfies the diffusion equation using the boundary condition of constant flux 
at the evaporation surface, that expresses as (reference[5] p. 61): —{C\D\ + 
(1 — C\)D2) '^r \ z =h= Ve, where C\ is the mean concentration of water in the 
layer defined as: C\(t m ) — j- J Q Ci{z,t m )dz at time t m . 

As a result an estimation of each quantity of equation [6] from the experi- 
mental results is required. Ci(t m ) is precisely estimated from the inset in figure 
[3J In the first drying regime, Ve can be measured from the mass variation 
with time in figure [5] ( Vg ~ 4 x 10~ 8 m/s). The values of the intrinsic coef- 
ficients are obtained from reference [TU]: ~ 50. The porosity K is roughly 
estimated to 0.40 for the random close-packed phase while f3 = 1/2 (reference 
[2"8]). We assume that the permeability of the gel keeps a constant value dur- 
ing the crack formation; this quantity is estimated using the Carman-Kozeny 
relation (k = 4.9 x 10~ 19 m 2 )[13 . Finally the elastic moduli of gels composed 
with various glycerol contents have been estimated from experimental results 
shown in figure [5jb) and (c) extrapolating values a short time after cracking. 
Consequently, the concentration distribution C\ and the pressure, P, are plot- 
ted along liquid pores at time t m in figure [5] for different glycerol contents. The 
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Figure 6: Concentration distribution in water, C±, (left) and normalized pressure 
distribution, PjP mi (right) along the pore liquid, for several initial glycerol 
contents, k. These quantities are plotted at time t m when the highest pressure 
value reaches P m at the evaporation surface. Inset: sketch of the pressure 
gradient in the liquid pore. 



flattening of the pressure gradient when the concentration of glycerol increases 
is clearly observable and is mainly attributed to the diffusion mechanism. One 
should note that the reduction in the capillary forces near the surface should be 
of minor importance as the surface tension from one mixture to another varies 
only from ~ 68mN/m to ~ 72mN/m [22J . These results highlight the effect of 
the addition of a non-volatile solvent to the pressure in the liquid pore since 
diffusion can extract liquid from the interior and flatten the pressure gradient 
in the pore liquid. The pressure gradient causes a shrinkage that is transmit- 
ted in the hard particles network. In the present paper, the adhesion with the 
substrate is assumed to be a constant and to play a determinant role as re- 
sponsible for the development of high stresses in the gel [5J [5] • Nevertheless the 
adhesion is seen as a necessary but not sufficient condition for cracking. Thus 
we eliminate cracking by modifying the pressure gradient and keeping constant 
the adhesion. The increase of glycerol content, through diffusion process, will 
imply that the resulting drying stresses will be distributed more uniformly in 
the gel. Consequently, the stress concentration in the gel network is limited and 
the crack formation can be inhibited. 

The effect of addition of glycerol content results in the flattening of the 



11 



pressure gradient in the pore liquid due to a diffusion process as depicted by 
equation H3 for a constant value of the elastic modulus. Moreover, the effect of 
addition of glycerol content results in the tendency of a decrease of the elastic 
modulus of the gel. As a consequence a decrease of the elastic modulus in 
equation [6] leads to a flattening of the pressure gradient, and, in that sense, 
emphasizes the effect of diffusion. 

5 Conclusion 

The incorporation of a small percentage of non- volatile cosolvent modifies the 
way the drying stresses relaxe in a colloidal gel. Experimental data show that the 
drying rate and cracking time is not significantly modified by incorporation of 
the cosolvent. However the presence of non- volatile decreases the stiffness of the 
gel at a given consolidation time and inhibits crack formation: incorporation of 
a concentration of glycerol greater than 10% in the colloidal dispersion results in 
a crack free drying gel. The mechanism is shown to be due to the combination of 
both flow driven by the pressure gradient and diffusion of solvent into a another; 
in this scenario the diffusion induced by evaporation causes a demixing effect of 
solvents. This results in a flattening of the pressure gradient in the pore liquid 
and consequently the distribution of the drying stresses is more uniform, limiting 
stress concentration in the gel network possibly cause of cracking. Moreover, the 
effect of a cosolvent added to a silica dispersion on the crack pattern has been 
highlighted in the case of glucose. As in the case of glycerol, the presence of 
glucose in the gel film reveals an increase of the crack spacing. However, contrary 
to the case of glycerol, cracks are not completely inhibited above a threshold 
glucose content: the crack spacing tends to a maximum value (typically 0.45cm 
in the same experimental conditions as for glycerol cosolvent). 

Finally, the effect of the flattening of the pressure gradient is surely not 
the only process leading to the diminution of the number of cracks in the gel 
film. Particularly chemical effects between glycerol and silanol groups covering 
the surfaces of silica particles could change the strength of the gel film. This 
process is suggested by aggregation rate measurements of silica particles with 
polvols[T51 and would be of interest for further works applied to gel film cracks 
induced by desiccation with chemical additives. 
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